function tbrillans_d(nu,b)

    use iasing_osse
    implicit none
    ! Arguments 
    real(kind=JPR8) :: tbrillans_d
    real(kind=JPR8), intent(in) :: nu,b

    ! Function variables
    real(kind=JPR8) :: nu0
    real(kind=JPR8) :: scale
 
    scale       = 100.0d+00
    nu0         = nu*scale
    tbrillans_d = Cst_hck*nu0/log((scale*Cst_2hc2*nu0**3/b)+1)
 
    return
end function tbrillans_d



